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ABSTRACT 

A theoretical explanation of some experimentally observed phenomena associated with the so-called Klebanoff 
modes is obtained by analyzing the flow over a finite thickness flat plate resulting from a small-amplitude distortion 
imposed on the upstream mean flow. The analysis shows (among other things) how the stretching of the vortex lines 
around the plate leads to streamwise vorticity at the plate surface, which then produces a streamwise velocity pertur- 
bation within the boundary layer that can be related to the experimentally observed Klebanoff mode. The complete 
evolution of this flow must be found by solving the boundary-region equations of Kemp (1951) and Davis and 
Rubin (1980), but a limiting analytical solution can also be obtained. Since the initial growth of the boundary-layer 
disturbance is nearly algebraic, our results demonstrate how the algebraically growing disturbances promoted by 
Landahl and others can be generated by a realistic external-disturbance environment. The relationship between these 
results and various bypass transition mechanisms is discussed. 


1. INTRODUCTION 

There are a number of different environmental disturbances that can cause transition to turbulence in boundary 
layers — with acoustic disturbances (Goldstein, 1985) probably being the most potent if they happen to be present in 
the flow. Another type of environmental disturbance that has received considerable attention in the literature is free- 
stream turbulence, which was actually studied by Dryden (1936), but the later investigations of Klebanoff (1971) 
seem to be more well known. 

Klebanoff measured the velocity fluctuations in a laminar flat-plate boundary layer due to externally imposed 
turbulence in a free stream (Reshotko, 1994) with nominally uniform upstream velocity. Figure lisa plot made by 
Klebanoff that has been shown around quite a bit, but, to our knowledge, has not actually been published in a journal 
article. He found, among other things, that the streamwise velocity fluctuations U M u' peak at an T| of about 2.3 and 
attain a maximum value of about one percent of the undisturbed free-stream velocity U M (wh ere t| is th e usual 
Blasius variable, i.e., the normal coordinate n divided by the boundary-layer thickness 8 0 = ,^/vx* /U M with x* de- 
noting the distance along the plate, and v denoting the kinematic viscosity). However, Klebanoff’ s most important 
observation was that the measured hot-wire signal passed through a low-pass filter at 12 Hz was almost identical in 
magnitude to the signal measured over all frequencies, indicating that most of the energy was at frequencies below 
12 Hz (Reshotko, 1994). This result has now been reproduced many times by Leventhal and Reshotko (1981), 
Kendall (1985 and 1990), Amal and Juillen (1978), Westin et al. (1994), and others. Kendall’s lateral correlation 
studies show these structures to be very narrow in the spanwise direction — of the order of twice the displacement 
thickness 5* in his particular case (where 8* ~ 1.7 8 0 ), while other experimenters find it to be nearly equal to 8*. 
Kendall called these low-frequency structures, Klebanoff modes, and that name seems to have stuck. 


NASA/TM— 1998-206298 


1 



Klebanoff (1971) proposed that these modes could be interpreted as a relatively small amplitude thickening and 
thinning of the Blasius boundary layer. Since perturbing the boundary-layer thickness S Q in the Blasius solution 
u = F b '(t|); r| ee n/8 0 (where F B denotes the Blasius function, a prime denotes differentiation with respect to the argu- 
ment, and all velocities are assumed to be normalized by UJ by a small amount, say 8 1# and expanding in a Taylor 
series gives 

« = F^(«/(5 0 + 5 1 ))=^(7|)-|-tjF|(i})+. . . . (1) 

This means that the difference between the actual streamwise velocity u and the Blasius solution should be propor- 
tional to tiF b , which leads to the excellent agreement shown in figure 2 — which is just a replot of the previous data. 

This idea was presumably based on an earlier proposal by Bradshaw (1965) and put on a more or less rigorous 
analytical basis by Crow (1966) to explain the fairly large thickening and thinning of the boundary layer observed 
by Klebanoff and Tidstrom (1959) in their famous peak- valley splitting work, under presumably disturbance-free 
conditions where the mean flow would be expected to be two-dimensional. Bradshaw (1965), who observed similar 
behavior at the National Physical Laboratory, proposed that the observed thickening and thinning could be produced 
by a small spanwise-periodic variation in the upstream mean-flow direction. 

Crow (1966) analyzed this phenomenon by considering a small-amplitude linearized sinusoidal perturbation of 
the flow over an infinitely thin flat plate, i.e., he considered the spanwise velocity perturbation (see fig. 3) 

w M = e cos 2 % z (2) 

imposed on a uniform flow over an infinitely thin plate, where e«l is a measure of the disturbance amplitude, X 
is its characteristic spanwise length scale, and z = z*/X. 

He assumed that the disturbance Reynolds number a = R^e (where R^ = U^X/v) is 0(1) and showed that the 
resulting streamwise boundary-layer velocity was of the form 


u 


= F B - — (ex sin 2 n z)77F B 


(3) 


where x = x*/X is the scaled streamwise coordinate, which means that the corresponding boundary-layer thickness 
would be 


5 = 5o(l + — xsinlnz +. . . . (4) 

V 2 J 


2. AN ALTERNATIVE PROPOSAL 

Crow’s solution leads to a steady streamwise vorticity disturbance at the surface of the plate. On the other hand, 
the most general vortical distortion that can be imposed on the flow consists of all Fourier components of an arbi- 
trary (unsteady) convected gust solution (Goldstein, 1978) of the linearized inviscid equations of motion. The 
streamwise vorticity component of the upstream motion was considered by Goldstein and Leib (1993a). However, 
here we look at another component of the upstream motion that might be visualized as the periodic steady wake flow 
shown in figure 4, and can best be described as a spanwise variation in the streamwise velocity, so that the upstream 
vorticity is initially normal to the plate. 

Our reason for considering this component is based on the following considerations. Firstly, the inviscid rapid 
distortion theory described below (see Hunt and Carruthers, 1990 and Goldstein, 1978, for additional details) shows 
that the inviscid vortex stretching around the plate tends to amplify the normal component of vorticity more than the 
other components. It also shows that, while the surface velocity remains bounded when the frequency of the distur- 
bance is 0(U oo /X), the zero-frequency (i.e., quasi-steady) component becomes logarithmically infinite at the surface 
of the plate— which means that the vortex stretching acts like a low-pass filter that preselects the low-frequency 
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component of the motion. Secondly, the solution to the linearized unsteady boundary-layer equations (Guliaev et al., 
1989) shows that the boundary layer on the surface of the plate also selectively amplifies the low-frequency portion 
of the spectrum with the small spanwise wavelength components undergoing the greatest amplification. The un- 
steady part of the flow actually moves out of the boundary layer at large values of x, (first into a kind of edge layer 
and then into the free stream) so that the boundary layer also acts as a low-pass filter for the free-stream turbulence. 

These various mechanisms are responsible for producing the low-frequency elongated structures in the bound- 
ary layer which we have referred to as Klebanoff modes. Since the external turbulent flow behaves more or less 
linearly to lowest order of approximation, we can probably gain some useful insight into the dominant motion in the 
boundary layer by considering the most amplified Fourier component of that flow, which leads to the upstream 
boundary conditions described above. 

While the inviscid vortex stretching around the plate tends to preselect the zero-frequency (i.e., co « U M /X,) 
component of the free-stream distortion, the boundary-layer flow itself actually selects the low-frequency compo- 
nents corresponding to cox*/U oo = 0(1), where co denotes the frequency, so that, while the external flow is quasi- 
steady, the boundary-layer motion is actually governed by the unsteady boundary-layer equations (Guliaev et al., 
1989). However, the purpose of this paper is to explain some of the relevant physical phenomena, rather than nu- 
merically reproduce the experimental results, and this can best be done by concentrating on the nonlinear phenom- 
ena and ignoring the unsteady effects. 

Moreover, Watmuff (1997) found that he was able to reduce the peak RMS velocity in the Klebanoff modes by 
about 50 percent and make his measured contours of constant RMS velocity much more uniform by reordering his 
screens based on laser scans of the uniformity of the mesh. He was also able to directly measure weak spanwise 
nonuniformities in the mean flow velocity upstream of the plate which had an elongated streakiness in the transverse 
direction. Watmuff concludes from this that, “The evidence is almost overwhelming that the nonuniformities in the 
free stream (and hence the Klebanoff modes in the boundary layer) originate from small nonuniformities in the 
porosity of the wind tunnel screens.” The wake-like distortion described above clearly provides an appropriate math- 
ematical representation for this type of nonuniformity. 

The stretching of the vortex lines around the plate not only amplifies the low-frequency component of the vor- 
ticity, but also produces the required streamwise vorticity at the plate surface as shown in figure 4. This streamwise 
vorticity then leads to a perturbation to the streamwise boundary layer velocity that is similar to the result obtained 
by Crow. Both these velocity perturbations increase more or less linearly with the streamwise distance x*, and the 
solutions are therefore not uniformly valid in the streamwise direction, i.e., the perturbation velocity will become of 
the same order as the base flow when ex = 0(1) (see (3)). Moreover, the continued downstream growth of the mean 
boundary layer eventually causes the boundary-layer thickness to become of the same order as the spanwise wave 
length of the upstream distortion. 

These linearized boundary-layer analyses can therefore only predict the initial stages of the boundary-layer 
development, and certainly not the inherently nonlinear processes that occur further downstream which may, as we 
shall see, be important for predicting the ultimate transition to turbulence. Moreover, the experiments show (as noted 
above) that the spanwise wave length of the dominant motion in the boundary layer (i.e., the Klebanoff mode) is of 
the same order as the boundary-layer thickness in the region where the measurements are carried out. 

3. ASYMPTOTIC STRUCTURE OF FLOW 

This leads us to consider the flow configuration depicted in figure 5. As in the Crow (1966) analysis, we allow 
the disturbance Reynolds number a to be 0(1), but now suppose that there is a wake-type velocity distortion im- 
posed on the flow at a distance 0 (A/e) upstream of the leading edge, where, as before, e is the amplitude of the 
distortion and X is its characteristic spatial scale. We also suppose that the plate thickness is 0(A.). Then in the as- 
ymptotic limit e — > 0, R^ — > °°, with o = eR^ held fixed, the flow divides itself into four distinct asymptotic regions. 
The first (Region CD in fig. 5) is a primarily inviscid region of dimensions 0(A,) surrounding the leading edge where 
the motion is a linear perturbation about a steady potential flow — which is the fundamental assumption for a branch 
of fluid mechanics now commonly referred to as linear-rapid-distortion theory (Hunt and Carruthers, 1990). Classi- 
cal rapid distortion theory (RDT) assumes that the disturbance length scale, X is much smaller than the mean-flow 
length scale. Here we follow the modem approach and allow these length scales to be of the same order. This theory 
has its origins in a marvelous paper by Sir James that he called simply “Drift” and published in the first issue of the 
Journal of Fluid Mechanics (1956). Lighthill considered the flow around an infinite circular cylinder and assumed 
that the imposed upstream distortion varied linearly with z, rather than being periodic as in the present study. 
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Below Region CD (i.e., the RDT region) is a region (Region © in fig. 5) where the flow is governed by the lin- 
earized laminar boundary-layer equations — as it is in the Crow (1966) analysis — with the resulting solution being 
somewhat similar to the one found by Crow — at least to the lowest order of approximation. As in Crow’s analysis, 
the linear boundary-layer solution breaks down when ex becomes large, and a new solution has to be obtained when 
x = ex is order one, i.e., in Region ®. The flow in Region © is now fully three-dimensional, even though the 
crossflow velocity w is still small, i.e., 0(e). This is because the streamwise length scale is now large enough to al- 
low a balance between the streamwise and spanwise convection terms within the governing equations. The flow in 
this region is then governed by what were originally referred to as the “boundary-region equations” by Kemp 
(1951), who was the first to use them, but that terminology seems to have gotten lost over the years. In any event, 
they are simply the Navier-Stokes equations with streamwise derivatives neglected in the viscous and pressure gra- 
dient terms and, more importantly, correspond to a rational asymptotic limit of the Navier-Stokes equations for the 
flows of the type being considered here. In the present paper, we resurrect the original terminology and refer to this 
approximate form of the Navier-Stokes equations as the boundary-region equations. 

Finally, the flow in the large outer Region © is very simple and corresponds to a slow viscous decay of the im- 
posed upstream velocity perturbation 


u = 1 + £u oo (z)e _x/0r + ... (5) 

over the long streamwise length scale which is balanced by a weak crossflow velocity 


(7 dz 

needed to satisfy continuity. As in the Crow analysis, we have put x = x*/X, y = y*/A,, z = z */X and have assumed, 
for simplicity, that u M = a cos 2k z, where a = constant, i.e., that the upstream distortion is purely harmonic. The 
result can easily be generalized to an arbitrary periodic function of z, but the final result is not as transparent, and the 
physical phenomena can probably be better illustrated by considering only a single harmonic component. 

The upstream boundary condition has to be imposed at a finite distance upstream of the leading edge, say L, on 
the long streamwise length scale x, but this corresponds to upstream infinity on the rapid inner scale x on which the 
potential flow field about the plate actually varies. So the upstream boundary condition is more or less unaffected by 
the presence of the plate in the strict asymptotic sense (but see Section 2.5 of Goldstein, 1978). Finally, the solutions 
in each of the individual Regions © to © must match together in the matched asymptotic expansion sense 
(Van Dyke, 1975). 


4. THE RAPID DISTORTION THEORY SOLUTION 
The flow in Region © (the RDT region) is given by the simple formula (Goldstein, 1978) 

u = VO + e[V(f) + au M (z)VA] (7) 

where 0(x,y) is the known potential for the mean flow about the plate (including the mean boundary-layer displace- 
ment effects, if necessary), a = e -L/a is a constant, u oo (z) is still assumed to be purely harmonic, and A(x,y) is the 
famous Lighthill (1956)-Darwin (1953) drift function to which the title of Lighthill’s (1956) paper refers. Its differ- 
ence between any two points on a streamline is equal to the time it takes a fluid particle to traverse that distance. It 
can be expressed as a simple integral (Lighthill, 1956; Goldstein, 1978; Darwin, 1953). Finally, the “perturbation 
potential” <j)(x,y,z) (which is produced by the interaction of the upstream distortion with the flow field of the finite 
thickness plate) is determined by the simple Poisson’s equation 

V 2 0 = -au DO V 2 A (8) 

subject to the boundary conditions 

(j) — > 0 as x — ^ -oo (9) 
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( 10 ) 


— + au^ — = 0 on surface of plate. 

<3n dn 

Its solution shows that the crossflow velocity w exhibits the well-known logarithmic singularity 

w->^-&i'F, (11) 

a 

as the mean-flow stream function T^x.y) goes to 0 at the surface of the plate (Lighthill, 1956; Goldstein, 1978) 
where a is a constant related to mean potential flow in the vicinity of the forward stagnation point. Since a similar 
analysis of the unsteady flow shows that the OfU^/X,) frequency components of the motion remain bounded at the 
surface of the plate, this result justifies our observations (in Section 2 above) about the selective amplification of the 
low-frequency component of the motion. 

The solution to the boundary- value problem (8) to (10) must, in general, be found numerically. However, we 
are mainly interested in the flow in the nonlinear Region ®, and the upstream boundary conditions for this region 
can be obtained from the asymptotic solution to the RDT problem as x — > °®. Fortunately, this result can be obtained 
independently of the upstream solutions (Goldstein and Durbin, 1980; Goldstein, Leib, and Cowley, 1992) by (ana- 
lytically) solving a simple equation with constant coefficients. The resulting solution, which is the same as in 
Goldstein, Leib, and Cowley (1992), depends on the upstream history effects only through the asymptotic drift func- 
tion A + + x, where A + is a suitable integral over the particle paths from upstream infinity to downstream infinity. The 
latter exhibits the logarithmic singularity 


A, — » A 0 (n)- — inn as n— >0 
a 

at the surface of the plate corresponding to the log singularity in (1 1). 


( 12 ) 


5. THE LINEAR BOUNDARY-LAYER SOLUTION 

Region © coincides with the boundary layer on the surface of the flat plate. The relatively mild singularity in 
the inner limit of the Region (D solution (11), which forms the outer limit of the Region © solution, can be smoothed 
out by viscous effects (Goldstein, Leib, and Cowley, 1992; Toomre, 1960) in Region ©. In fact, the mean boundary 
layer is again of the Blasius type sufficiently far downstream in the flow, and an analysis similar to that of Crow 
(1996) shows that the distorted streamwise velocity is now given by 

u Fg + + 0(£x) (13) 

to lowest order of approximation. This result differs from Crow’ s because the gauge function ex €n (ex) is larger 
than his by the factor in (ex) in order to accommodate the log singularity in the external flow. However, the next 
order (i.e., 0(e)) terms in (13) would be significantly different from those obtained by Crow (1966) (whose lowest 
order perturbation is 0(ex)). 


6. THE LARGE-DISTURB ANCE-REYNOLDS -NUMBER LIMIT OF THE BOUNDARY-REGION 
EQUATION SOLUTIONS AND THE FORMATION OF TURBULENT SPOTS 

The asymptotic solution (13) provides the upstream matching condition (i.e., the upstream boundary condition) 
for the flow in the larger region downstream (Region © in fig. 5). While the solution in this region is parabolic and 
can therefore be found by downstream marching, it still constitutes a relatively difficult numerical problem (Davis 
and Rubin, 1980; Wundrow, 1996). It therefore seems desirable to obtain some simplified analytical results. This 
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can be done by considering the large- and small-disturbance-Reynolds-number limits. In the small-disturbance- 
Reynolds-number limit, the boundary-region equations become linear. Here we discuss only the large-disturbance- 
Reynolds-number limit in which g = eR^ goes to infinity (as was done in Goldstein, Leib, and Cowley, 1992, and 
Goldstein and Leib, 1993a, b). 


In the large-G limit, the downstream nonlinear flow splits into two layers when Xj = — fn 0(1) where 

the new gauge function S(a) is determined by i = ^Jcrfn^A/Sj : a viscous wall layer in which the flow is governed 
by the three-dimensional, boundary-layer equations with no pressure-gradient terms — rather than by the boundary- 
region equations and an outer inviscid vorticity layer induced by the log singularity in the upstream RDT solution. 
The scaled crossflow velocity W = awMn(A/5) in this latter region is determined by the inviscid Burgers equation 
(see fig. 6 ) 


dW <9W 

— + W^— = 0, 

dz 


(14) 


which can, of course, be solved analytically. This solution can then be used to obtain the appropriate outer-edge 
boundary condition for the viscous wall layer whose thickness <5 is now smaller than X by a factor of the square root 
of g (the large disturbance Reynolds number) times €nXJ S (Goldstein, Leib, and Cowley, 1992). The boundary 
condition is that the streamwise velocity U goes to one, and the scaled crossflow velocity W goes to 

W— >f(z-x 1 W) as y = n/<5->°o (15) 


where 


f(0 = -uL(0. 

u M z is the imposed upstream velocity perturbation discussed above, and the primes denote differentiation with re- 
spect to the total argument (Goldstein, Leib, and Cowley, 1992). 

The solution to this boundary-value problem can be expressed in terms of the Blasius solution by 

U = F£(77) (16) 


where 77 is now given by 


77 = 


(l + xif')y 


3f' 


(1 + Slff-l 


(17) 


rather than by the usual Blasius variable. Timoshin (1992) recently used a similar transform to solve a very different 
boundary-layer problem. Since rf now depends on z, U also exhibits this dependence, and the solution is therefore 
directly interpretable as a localized thickening and thinning of the Blasius boundary layer. It is a natural extension of 
the linearized Crow -type solution into the nonlinear region. In fact, the large-G limit of the linear boundary-layer 
solution (13) can be recovered by expanding F B and 77 in a Taylor series for small X] — so the solution (16) clearly 
matches onto the appropriate linear solution in Region © as — » 0. 

Of course, the solution to Burgers equation (14) eventually develops a singularity at a finite downstream posi- 
tion and is discontinuous or multivalued downstream of that point. This introduces a line of singularities into the 
vorticity-layer solution and a surface of discontinuity downstream of that line (Goldstein, Leib, and Cowley, 1992). 
(See fig. 7.) A new inviscid solution, which accounts for streamline curvature effects, can then be constructed to 
eliminate the line singularity, and (downstream of that) an inviscid solution that brings in pressure-gradient effects 
can be constructed to eliminate the surface discontinuity (Goldstein and Leib, 1993b). This is all rather intricate, but 
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the important thing to note is that this solution is still inviscid and has the general structure shown in figure 7 
(Goldstein, Leib, and Cowley, 1992; Goldstein and Leib, 1993b). However, the downstream pressure-gradient solu- 
tion imposes its strong pressure gradients on the underlying boundary layer and causes it to develop a singularity at 
a finite downstream position. This singularity is an indication that the boundary-layer approximation breaks down, 
and a new (less approximate) set of equations has to be used to describe the flow. Since the spanwise length scale 
decreases and the boundary-layer thickness increases in the vicinity of the singularity, this probably turns out to be 
the full boundary-region equations. 

The boundary-layer velocity profile eventually becomes much flater and can even become inflectional in the 
viscinity of the singularity. (See Goldstein, Leib, and Cowley, 1992.) The flow therefore becomes susceptible to 
rapidly growing Tollmien-Schlichting waves and possibly even Rayleigh instabilities which could ultimately lead to 
the formation of turbulent spots. It therefore provides a possible “bypass” mechanism that could lead to transition 
without involving a more global Tollmien-Schlichting wave system (but see below). 

Figure 8 is a plot (taken from Goldstein, Leib, and Cowley, 1992) of constant streamwise velocity contours in 
the crossflow plane calculated at various streamwise locations with u M set equal to cos z. The initial spanwise varia- 
tion in the linear regime is, of course, sinusoidal, and the subsequent effect of the nonlinearity is to produce a sharp 
focusing along the lines of maximal thickness, with the attendant introduction of progressively shorter streamwise 
length scales into the flow. 

Since the large streamwise velocity perturbations are concentrated in narrow spanwise regions, the spanwise 
average RMS velocity 


J. 

it 1 2 

iJ[U(x,,Iu)-F4(7,)] 2 dz (IB) 

0 J 

(which roughly corresponds to the local RMS velocity measurement in the real flow) will be much smaller than the 
maximum velocity fluctuation. This is shown in figure 9, in which the solid curve is the peak value (relative to T)) of 
the RMS velocity (18), while the dashed curve is the peak (with respect to Tj) of the maximum velocity fluctuation 



|u(x 1 ,77,0)-F^(?7)| 

with U determined by (16) and (17) in both cases. Notice that while the RMS velocity fluctuation is only a little over 
10 percent at xi — 0.75, the maximum velocity fluctuation is nearly equal to 50 percent at this point. Figure 8 shows 
that this is large enough to produce localized nonlinear behavior, while the straight-line behavior of the solid curve 
suggests that the RMS velocity continues to behave linearly. 

The experiments (see comprehensive summary in table III of Westin, et al., 1994) show that the Klebanoff 
mode RMS velocity fluctuations can easily exceed 10 percent. The present results, therefore, suggest that the corre- 
sponding maximum velocity fluctuations will then be large enough to produce the localized breakdown of the 
boundary-layer solution alluded to above, even though the RMS velocity fluctuations are still relatively small. Need- 
less to say, nonlinear theory will have to be used to calculate the flow. 

Figure 10 is a flow visualization by Kendall (1985) using a smoke-wire technique to view the crossflow plane 
with a camera located 3m downstream of the test section. The smoke wire was located at a Blasius t| value of about 
1.5. This figure suggests that the crossflow distortion is focused in rather narrow spanwise regions in agreement with 
the theoretical predictions. 

7. APPLICATION OF SOLUTION TO KENDALL’S EXPERIMENT 

Figure 1 1 shows some data taken by Kendall (1985, 1990) in an experiment similar to Klebanoff s. The closed 
symbols show the agreement of the mean velocity distribution with the Blasius profile. The open symbols are the 
filtered signal at the Tollmien-Schlichting wave frequency. This disturbance is clearly small (note the 20-fold mag- 
nification) and concentrated at the outer edge of the boundary layer — definitely not a T-S mode. The remaining 
curve represents the broad-band signal corresponding to the low-frequency Klebanoff mode. The solid curve is the 
rjFg mode shape proposed by Klebanoff and put on a more rigorous analytical basis by Crow (1966). Kendall’s 
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turbulence levels are higher than those of Klebanoff by a factor of about four or five, but the results are the same; the 
agreement with Crow’s mode shape is excellent, except perhaps at the outer edge of the boundary layer where the 
data approaches a finite value while Crow’s analytical distribution goes to zero. Figure 12 is a replot of the data in 
Figure 11, but the solid curve is now obtained from (18) with U determined from (16) and (17). The agreement is 
now slightly better than it was with Crow’s linear mode shape. 

While Kendall did not observe any Tollmien-Schlichting waves in the region where the Klebanoff modes first 
become predominant, he did detect T-S wave packets further downstream. These wave packets frequently exhibited 
nonlinear amplitude dependence on the free-stream turbulence level, but much more interestingly, appeared to be 
highly intermittent — even though Kendall’s measurements showed the free-stream turbulence to be relatively time- 
stationary (Kendall, 1985). 

A possible explanation for this, as well as for the nonlinear dependence on free-stream amplitude, is that the 
upstream Klebanoff modes act as receptivity sites for the free-stream turbulence (Goldstein, Leib, and Wundrow, 
1992). At the low Mach numbers of Kendall’s experiment, there is probably a significant mismatch between the T-S 
wave length and the wave length of the free-stream turbulence — with the latter being considerably longer than the 
former. It is therefore necessary that the boundary-layer flow contain some sort of short streamwise length scale that 
can “scatter” the long wavelength free-stream disturbance into the shorter T-S waves. It has been shown (Goldstein, 
1983; Goldstein et al., 1983; Goldstein and Hultgren, 1989; Goldstein, 1985) that (1) imperfections in the plate sur- 
face, (2) the sudden change in the curvature at the junction between the leading-edge region and the flat portion of 
the plate, and (3) the rapid growth of the laminar boundary layer at the leading edge can all produce the requisite 
short streamwise length scale. Our proposal here is that the Klebanoff modes can also fulfill this role with the short 
streamwise length scales resulting from the previously shown focusing effect that occurs in the large-o asymptotic 
solution. 

Since the Klebanoff modes are, in reality, not very steady, but rather meander across the plate at low frequen- 
cies, i.e., they are governed by the unsteady boundary-layer equations, the more or less localized receptivity centers 
would also be expected to move around and thereby cause the T.S. wave packets to appear to be intermittent to a 
downstream observer with a fixed probe. 


8. NUMERICAL SOLUTIONS FOR THE BOUNDARY-REGION EQUATIONS AND ALGEBRAIC 

DISTURBANCE GROWTH 

The appearance of logarithmic gauge functions in the large o asymptotic expansion for the solution in Region ® 
(see Goldstein, Leib, and Cowley, 1992) suggests that this expansion is not very robust. The only recourse is to 
solve the full (or perhaps the linearized) boundary-region equations numerically. Here we discuss the full nonlinear 
case. The relevant initial conditions are determined by requiring that the solutions match onto the asymptotic solu- 
tions in the upstream RDT/linear boundary-layer region (Region (D and © in Fig. 5). Since, as we have seen, the 
streamwise boundary-layer velocity grows nearly linearly with x (actually it grows like x Tn x), the initial velocity in 
Region © will have to exhibit the same nearly algebraic growth. (This is, of course, also true for the large-o asymp- 
totic solution described above, which actually exhibits purely linear growth to lowest order of approximation.) 

This linear or nearly linear growth corresponds to the algebraic growth mechanism originally discussed for 
boundary layers by Hultgren and Gustavsson (1981), subsequently promoted by Landahl and his students (Landahl, 
1990; Russell and Landahl, 1984) as a possible explanation for certain types of bypass transition, and now appears 
to be gaining increasing acceptance by the transition community as a viable alternative to instability wave mecha- 
nisms. We have applied it to spatially growing disturbances (in Goldstein, Leib, and Cowley, 1992, and Goldstein, 
Leib, and Wundrow, 1992) and, more importantly, shown how these latter disturbances can be generated by realistic 
disturbances in the environment (i.e., we have shown in Morkovin’s words that these algebraically growing distur- 
bances are “environmentally realizable”). Our results also suggest that it is probably necessary to account for non- 
parallel flow effects when dealing with disturbances that are both environmentally realizable and algebraically 
growing. 

Figure 13 shows some preliminary solutions to the boundary-region equations taken from Wundrow (1996). 
These results roughly match onto the upstream RDT solution with u M still given by cos z. However, the initial con- 
dition has now been adjusted to make the velocity profiles more inflexional. The figure again shows contours of 
constant streamwise velocity in the crossflow plane at various streamwise locations. As expected, the initial flow is 
relatively two-dimensional with the perturbation varying sinusoidally in the z-direction. The figure also shows a 
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fairly rapid departure from the initial two-dimensional state that becomes more nonharmonic with increasing down- 
stream distance. The relative distortion becomes maximal somewhere downstream of t, = 1.3 (i.e., x = 0.3a/e) and 
then rapidly decays out with further increase in downstream distance due to the action of viscosity. So, unlike the 
large-disturbance-Reynolds-number limit that we discussed above, the boundary layer does not develope a singular- 
ity in this case. 

The ideas discussed in this paper can be extended to the case where the upstream distortion is unsteady. The 
boundary layer will filter out most, but not all, of the unsteady components of the motion, and the flow in Region ® 
will be described by the unsteady boundary-region equations. The relevant solution to these equations requires ex- 
tensive numerical computation, which can, however, be greatly simplified by considering the small-disturbance- 
Reynolds-number limit. This has been done in Leib, Wundrow, and Goldstein (1997). 


9. EFFECT OF KLEB ANOFF MODES ON INSTABILITY WAVES 

At very low free-stream turbulence levels (~0.1%), even the Schubauer and Skramstad (1948) experiments 
show that the nonlinear breakdown of Tollmien-Schlichting waves are directly responsible for the appearance of 
turbulent spots. The situation is less clear at intermediate turbulence levels (of the order of 1 percent or so). Figure 
13 of Westin et al. (1994) shows that there is enhanced streamwise amplification of the boundary-layer fluctuations 
within the band of frequencies corresponding to unstable Tollmien-Schlichting wave growth. However, these fluc- 
tuations tend to be broad band, rather than being concentrated in a narrow range about the most unstable frequency 
as in the low-turbulence-level experiments. But follow-on experiments by Boiko et al. (1994), in which discrete 
frequency Tollmien-Schlichting waves were artificially introduced into the flow, show that sufficiently high levels 
of free-stream turbulence (and, therefore, presumably the Klebanoff modes) can produce significant transfer of en- 
ergy between frequencies within the unstable bands for Tollmien-Schlichting waves. It may therefore be possible 
that unstable Tollmien-Schlichting-type waves growing on a base flow computed from the unsteady boundary-re- 
gion equations will exhibit a sufficiently broad band- width to explain the experimental observations. But in order to 
verify this, it is necessary to investigate the stability of this unsteady base flow, which is currently being done by the 
present authors. However, even without doing the calculation, it is clear that certain portions of an initially 
undistorted T.S. wave would grow faster than others as it propagates downstream, since certain portions of the base 
flow would be more unstable than others. In fact, the rapidly growing localized Tomlmien-Schlichting waves and/or 
Rayleigh instabilities that we discussed in conjunction with the large turbulent Reynolds-number, Klebanoff-mode 
solution could be part of a more global T.S. wave system that impinges on the locally destabilized region from up- 
stream. This would certainly enhance the local growth of the T.S. wave over the two-dimensional case, even though 
the average growth could be much smaller than in the two-dimensional case. 

The various mechanisms discussed in this paper could be operative in different flow regimes up to free-stream 
turbulence levels of about 5 percent. Beyond that point, it is impossible to detect any turbulent spots or, for that mat- 
ter, any other discernible features of the flow. 


10. CONCLUDING REMARKS 

In all of the experiments of which the authors are aware, no matter what the free-stream turbulence level, the 
dominant streamwise velocity fluctuations are of the Klebanoff type, i.e., they are of very low frequency, have a 
peak value at a Blasius q -value of about 2.3, are very elongated in the streamwise direction, and are very narrow in 
the spanwise direction. We have noted that this is consistent with the predictions of inviscid rapid- distortion theory, 
which shows that the vortex stretching around the plate selectively amplifies the low-frequency component of the 
motion. We also pointed out that it is consistent with the solutions to the linearized unsteady boundary-layer 
equations which, for the convected (i.e., pressureless) free-stream disturbances that describe the weak free-stream 
turbulence, grow linearly with the streamwise coordinate in the low-frequency limit and move out of the boundary 
layer (first into a kind of edge layer and then into the free stream) in the high-frequency limit; so that only low-fre- 
quency disturbances remain in the boundary layer at sufficiently large distances downstream. 

Moreover, since it turns out that the streamwise velocity perturbations are driven by the spanwise gradient of 
the spanwise velocity, the streamwise velocity growth rate initially increases with decreasing spanwise length scale. 
This trend is eventually reversed by viscous effects in the downstream region where the boundary-layer thickness 
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becomes of the same order as the spanwise length scale, and, in the most general case, the motion is governed by the 
unsteady boundary-region equations. But even if this selection did not occur, the free-stream turbulence and/or the 
mean-flow distortion which drive the Klebanoff modes would eventually run out of energy when the spanwise 
length scale became too short. These two factors combine to form a selection mechanism for the spanwise length 
scale, with one or the other becoming dominant, depending on the nature of the free-stream turbulence. This might 
explain the continuing argument among experimentalists as to whether the spanwise length scale is set by the free- 
stream turbulence or by the boundary layer itself. 
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Tf| = n/5 0 r) = y/8 0 


Figure 1 . — Streamwise velocity fluctuations in laminar Figure 2. — Comparison of Klebanoff's data with Taylor 

boundary layer measured by Klebanoff (unpublished series expansion of Blasius solution. Same data as 
data) x| = n/8 0 where n is normal distance to wall, in Figure 1 . 

8q = /0c* v/U oo is the Blasius boundary layer thickness. 

The streamwise distance measured from leading 
edge x* is equal to 22 inches. 



Figure 3. — Flow configuration analyzed by Crow (1 996): A. = spanwise wave 
length of distortion; w^ = imposed crossflow velocity. 
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Figure 4. — Schematic view of stretching of vortex lines around plate: 
\ = spanwise length scale of distortion, e = characteristic scale of 
imposed streamwise velocity distortion, = imposed vorticity 
distortion, c*> = local vorticity lines. 



Figure 5. — Overall flow configuration showing various asymptotic regions 
considered in the analysis. 
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Figure 6. — Asymptotic flow structure in boundary layer in high disturbance 
Reynolds-number limit. 



Figure 7. — Singularity structure of the nonlinear inviscid flow from 
Goldstein, Leib, and Cowley (1992). 
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Figure 9. — Comparison of RMS and maximum velocity 
fluctuations in boundary layer as a function of the scaled 
streamwise coordinate x 1 . Solid line; peak value of 
u'( x 1s ti), dashed line; peak value of |U( x lP y , o) — F' b (t|)|. 



Figure 10. — Flow visualization of boundary layer beneath a turbulent stream (from Kendall, 1985). 
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Figure 1 1 . — Streamwise velocity fluctuations in boundary layer as measured by 
Kendall (1985). Turbulence is generated by a 6 cm grid /\/U«,x*/v = 717. 



0 2 4 6 8 1( 

Figure 12. — Comparison of RMS velocity profile (18) at x 1 = 0.27 with 
experimental results of Kendall (1985). 
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